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We argue that, due to the specific form of the interaction potential between electronic 
guiding centers, a bubble crystal (BC) with basis may be energetically more favorable 
than the usual bubble solid. This new BC has well-defined normal eigenmodes, and 
remains energetically favorable when the effects of finite sample thickness and of screening 
by electrons in lower LLs are taken into account. 
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1. Introduction 

Quantum phases of the two-dimensional electron gas in higher Landau levels have 
been the subject of intense study in recent years. Hartree-Fock studies 1 have shown 
that, for small partial filling factors v* = v — 2n w 0.1 — 0.2 (with v the total 
filling factor and n the Landau level index) , the electrons form a triangular Wigner 
crystal, while for v* close to 1/2 the ground state of the system is a stripe state. 
Between these two regions, a new crystalline solid with more than one electron per 
lattice site, termed bubble crystal, turns out to be energetically more favorable. 
In this paper, we investigate the interaction of electrons within individual bubbles, 
and argue that a bubble lattice with basis may lead to energetically more favorable 
ground states of the two-dimensional electron system, the properties of which we 
shall briefly describe. 

2. Wigner crystal and bubble phases in higher Landau levels 

We shall start from the expression of the Hartree-Fock energy of the partially filled 
nth Landau level, which is given by: 



where Vhf is the Hartree-Fock interaction, and the order parameter p(q) is 
the Fourier transform of the guiding center density, and is given by p(q) = 
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n(q) /[e~ q2 / 4 L n (q 2 1 2 /2)] . In this last equation, n(q) is the Fourier transform of 
the electron density, which we shall approximate by n(r) = J2i m l<£nm( r — Ri)| 2 , 
where y„ m (r) is the noninteracting wavcfunction of angular momentum to and Lan- 
dau level index n, and where the summation extends over all the bubbles located at 
the lattice sites R^ of a triangular Bravais lattice and all the electrons within each 
bubble. The Hartree-Fock interaction potential in Eq. (1) consists of the sum of a 
Hartree and Fock parts, which are given by: 

^(q)-^e-^ 2 / 2 [ L „ ((Z ^/2)l 2 , (2) 



eq 



(3) 



In the above expressions, e is the dielectric constant of the host semiconductor, and 

L n (x) is the nth Laguerre polynomial. 

It can easily be checked that the projected density p(q) can be written in 

the form p(q) = J2m Pm(<l) J2i e iq Ri i tne partial guiding center density 

p m (q) = n m (q)/[e-^ 2 /4 in((? 2 £ 2 /2)] (hcre ~ m(q) = j rfr |^ ro ( r )|2 e -iq-r is the 

Fourier transform of the density of electrons in state to) . Using the above expression 
of /?(q) into Eq. (1) allows us to rewrite the cohesive energy Ehf in the form: 

E ™ = \ E E - R ^) + E E (0) , (4) 

i^j m,m' i m<m' 

where we introduced the interaction potential U m ,m' between electrons in states m 
and m', which in real space is given by: 

U m , m '(r) = J Pm(q)V H F(q)p m '(ci)e^ r . (5) 

It can be verified that the above interaction potentials, when used in Eq. (4), give 
cohesive energies that are in excellent agreement with the Hartree-Fock calculations 
of Refs. 2 and 3. 



3. Lattice with basis: possible new bubble states 

A plot of the interaction potential ?7oi( r ) (through which the two electrons on a 
given bubble interact) reveals the existence of a local minimum at r = ro = 1A8£ 
(see Fig. 1), which might favor a finite separation between the guiding centers of the 
two electrons in a given bubble. Aiming at finding an energetically more favorable 
ground state of the bubble type for our two-dimensional electron gas, the simplest 
trial ground state one can think of consists of a triangular lattice of bubbles with two 
electrons per lattice site, but with the electrons guiding centers r m i at lattice site 
Rj separated from each other by a distance r , which may be obtained by writing 
r m i = Hi + (to — ^)rox (with to = 0, 1). Such a ground state considerably lowers 
the cohesive energy of the two-electron bubble crystal, as can be seen in Fig. 1. 
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Fig. 1. Left panel: Interaction potential U~oi (in units of e 2 /el) vs. r (in units of £). Right panel: 
Cohesive energies (in units of e /el) of the Wigner crystal (thin solid line), ordinary bubble crystal 
(dotted line), stipe phase (dashed line) and of the bubble crystal with a basis (thick line). The 
cohesive energy of the stripes is calculated as in Ref. 3. Both figures are for Landau level n = 2. 

4. Normal modes and stability of the bubble crystal 

We now want to check the stability of our bubble crystal with basis. We shall 
do so by finding its normal vibration modes, and verifying that they are all 
well defined. The derivation of these normal modes will proceed in a standard 
way as follows. To fix ideas, we shall consider the simplest case of the two- 
electron bubble crystal in the n = 2 Landau level, and write for the electron 
guiding centers r mi (to = 0, 1) at a given lattice site Ri the decomposition 
r mi = Hi + (to — ^)r x + u mi , where u mi is the displacement of the m-th 
electron from the equilibrium lattice position R^. Expanding the cohesive en- 
ergy E HF = \ E 1?ij E,», m - ^^'f 1 ".. - V,) + E, E m <m' ^.m' ((» ~ «!>«X + 

u m i — u m 'i) to second order in the small displacements \i m i leads to the elastic 
energy E el = \ £\ ^ J2 m , m ' u ma(Ri)$X ( R » ~ R-j)«m'/3(Rj) » where the elastic 
matrix $ is given by ^"^'(Ri - Rj) = 6ij6 mm > Y,k, m " 9 a dpU mm »(Kk + (m - 
m")rox) — d a dpU mm i{TLi — R^ + (m — m')r x) . In presence of a magnetic field, 
the equation of motion for the m-th electron at lattice site Ri has the form 
m*^u ma (R t ) = -E,>/*X'(Ri - R>m'/3(R-j) - ^ a/ sf t u m0 (-Ri) (m* is 
the effective mass of the electron in the host semiconductor and e a p is the com- 
plete antisymmetric 2D tensor). We shall seek a solution to the above equation 
of motion that represents a wave with angular frequency uo and wavevector q, 
u mce (Ri) = A ma {q) e *(q- R i- w *) . Substituting the above expression into the equation 
of motion, and solving the resulting secular equation leads to the normal eigenmodes 
plotted for v — 4.20 in Fig. 2 (left panel). As can be seen, the normal modes are all 
real, and consist as is well-known of two magnetophonon modes with characteristic 
dispersion u>(q — > 0) ~ q 3 ^ 2 , and two magnetoplasmon modes with uj(q — > 0) = uj c 
(uj c = eB/m*c being the cyclotron frequency). 

We have also investigated the combined effects of screening by electrons from 
lower LLs (encoded in a wavevector-dependent dielectric constant e(q)) and of a 
finite sample thickness on the cohesive energy of the bubble crystal with basis. The 
right panel on Fig. 2 shows the cohesive energies of the various phases vs. v* (at LL 
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q/(7c/a) v* 

Fig. 2. Left panel: Magnetophonon and magnctoplasmon dispersion curves along the boundary of 
the irreducible element of the first Brillouin zone of the BC with basis at v = 4.20 (a is the bubble 
lattice constant and uiq = e 2 /he£). The left hand side frequency scale is for magnetophonons, and 
the right one is for magnctoplasmons. Right panel: same as the right panel in Fig. 1, but with the 
g-dependent dielectric constant e(q) and with the finite sample thickness parameter A = 1. 

n = 2), with the Coulomb potential Vo(q) = (2ire 2 )/eq in Eqs. (2)-(3) replaced by 
v(q) = (2we 2 )e~ Xqe / e(q)q. The parameter A models a finite thickness sample, 4 and is 
generally taken to be of order unity. Here we use A — 1, while for e(q) we shall use the 
following expression, due to Aleiner and Glazman, 5 e(q) = e(l + -^[l — Jq (<7_R C )]^ , 

with as = h 2 e/m*e 2 the effective Bohr radius and R c = \j2n + 1 1. We find that, 
although screening and finite sample thickness shift the cohesive energies up, the 
relative boundaries between the various phases remain practically unchanged. 

5. Conclusion 

In conclusion, in this paper we have argued that a new bubble crystal with a basis 
may be an energetically more favorable bubble ground state of the 2D electron 
system. We have shown that this new bubble state with basis is stable against small 
fluctuations of the guiding center positions, and remains energetically favorable in 
finite thickness samples and if we take screening by lower LLs into account. More 
detailed investigations of other possible crystal structures for the Wigncr crystal 
and bubble phase (along with experimental consequences) are under way and will 
be discussed elsewhere. 
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